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Recent measurements of a peak in the angular power spectrum of the cosmic microwave back- 
ground appear to suggest that geometry of the universe is close to being flat. But if other accepted 
indicators of cosmological parameters are also correct then the best fit model is marginally closed, 
with the peak in the spectrum at larger scales than in a flat universe. Such observations can be 
reconciled with a flat universe if the fine structure constant had a lower value at earlier times, which 
would delay the recombination of electrons and protons and also act to suppress secondary oscil- 
lations as observed. We discuss evidence for a few percent increase in the fine structure constant 
between the time of recombination and the present. 

PACS Numbers : 98.80. Cq, 95.35-fd 



I. INTRODUCTION 

Cosmologists have for many years struggled to find a model of the universe consistent with all the available evidence. 
Recently, many observations have pointed to the universe being spatially flat with cold dark matter (CDM) making up 
approximately a third of the critical density, and the remainder dominated by a component with a negative equation of 
state such as a cosmological constant, A. This confluence of evidence, which includes measurements of the expansion 
and acceleration of the universe, bounds on the age of the universe, and constraints from large scale structure and 
galaxy clusters, has been called 'cosmic concordance' 0. 

The restriction to spatially flat models was originally motivated by theoretical arguments, in particular to be 
consistent with inflationary models of the early universe. However, measurements of the anisotropy in the cosmic 
microwave background (CMB) have given strong observational support to this assumption. The position of the first 
Doppler peak in the CMB power spectrum is sensitive to the spatial geometry and the epoch of last scattering of the 
CMB photons relative to the present age of the universe Recent measurements by the balloon borne detectors 
BOOMERanG [D (B98) and MAXIMA [Q (M99) indicate an increase in power at angular scales of approximately 
two degrees, very close to the expected position of the first Doppler peak in spatially flat models. 

While providing tentative conflrmation of a flat universe, the new data has also raised many questions as to the 
precise viabihty of the concordance [|| ||. In particular, the position of the flrst peak as detected by BOOMERanG 
(^poak = 197 ± 6) appears to be at slightly larger scales than expected in generic flat models — a conclusion which 
is only slightly weakened by the inclusion the less sensitive MAXIMA data |^. In addition, both data sets indicate 
secondary Doppler peaks that are much less pronounced than expected in models compatible with primordial Big 
Bang Nucleosynthesis (BBN) 0,|l|- There are, of course, numerous possible explanations for these discrepancies. 
For example, the shift of the peak to larger scales could indicate a universe which is slightly closed [^|J|], while the 
suppression of the second peak could be evidence that the baryon density is higher than has been indicated by BBN 
However, these solutions require either giving up the elegance of spatially flat models or run into direct conflict 
with other cosmological measurements, particularly those related to BBN. 

Another possible solution which addresses both unexpected features of the data is to delay the epoch of last 
scattering. This would increase the size of the sound horizon at last scattering and shift the flrst peak to larger scales 
while keeping a spatially flat universe. It would also simultaneously increase the density of baryons relative to that of 
the photons during the epoch of last scattering, suppressing the amplitude of the second peak. Within the standard 
framework, it is rather difficult to change the time of decoupling since it would require a mechanism, astrophysical p^ ] 
or otherwise p"3[-p^, which could delay the formation of neutral hydrogen. Peebles et al. ||l2| suggested that non-linear 
structures at extremely high redshift could act as a source for Lyman-a photons which photo-ionize the hydrogen. 
However, this is very unlikely in standard adiabatic models for structure formation, though it is a possibility if the 
initial fluctuations were non-Gaussian. 

In this paper, we focus on another possible mechanism for delaying photon decoupling: that the electrons and 
protons might have been more weakly bound at high redshifts than they are today. In particular, we consider whether 
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the observations contain evidence for a running of the fine structure constant, a = / {AttTlc) , between the time of 
recombination, where the CMB last scattered, and the present epoch. Changes in a modify the parameters governing 
recombination which, depending on the sign, can lead to early ((5a(troc) = Aa(ii.cc)/a > 0) or delayed 

{5a{t-[cc) < 0) recombination. Here, we have defined A(t) = a{t) — a(to), where t is cosmic time and we denote tg, 
tree and tnuc throughout as the times of the present day, recombination and nucleosynthesis respectively. 

Such variation of physical constants has been the subject of much attention, both observational and theoretical, 
in recent years. The theoretical motivation comes from String or M-theory in models where there are compact extra 
dimensions. These extra dimensions may have either stabilized before recombination ((5a(iroc) = 0) or they may still 
be rolling down their potential, causing all the coupHng constants to vary (i5a(trcc) 7^ 0). The standard way in which 
to do this is using a scalar field known as the dilaton, but no stabilizing potential has ever been derived from anything 
which could be described as a candidate fundamental theory; all proposed stabilizing mechanisms appear to be ad 
hoc. Slow variation in a could, therefore, be considered at some level as a prediction of fundamental theory. It has 



also been pointed out |18 19[| that models in which a varies can be thought of as models with a varying speed of 

light [|o|,|il. 

Of course, limits exist on changes in a due to various terrestrial, astrophysical and cosmological arguments. Terres- 
trial limits come primarily from elements which have long-lived (3 decay, atomic clocks and the OKLO natural nuclear 
reactor [|2^, for which the limit is —0.9 x 10~^ < 5a < 1.2 x 10""^ over a time period of around 1.8 billion yearsj^, 
although this is model and theory dependent since the limits are sensitive to possible simultaneous variations in other 
coupling constants. Cosmological limits come from the Helium abundance in BBN and quoted limits can be expressed 
roughly as |(5a(inuc)| < 10^^ — 10^^ [p3|-p5| at z w 10^ — 10^°, although this is again highly model dependent; we 
shall return this issue in a detailed discussion below. Astrophysical limits arise from systems which absorb quasar 
emissions over a wide redshift range of z ~ 0.1 — 3 p6t . After many years of deriving upper bounds on (5a, a statistical 
detection oi 5a — (— 1.1 ± 0.4) x 10~^ has been claimed recently due to measurements of relativistic fine structure 
in absorption systems in the range 0.6 < z < 1.6, with more data soon to be published. 

In what follows we will work within the framework of spatially flat (f^tot = 1) ACDM models with matter density (in 
units of the critical density pcrit = 3i?g /SttG) given by il^ and cosmological constant density fl/^ = 1 — ^m- Similarly, 
the baryon density is denoted fib and the rest of the matter is assumed to be dominated by CDM, with no hot dark 
matter (HDM) component, 17^ = 0. The Hubble constant will be parameterized by Hq = lOOft-kmsec"^ Mpc~^. We 
shall assume that the initial scalar fluctuations that are measured were created during an epoch of cosmic inflation, 
and that they are almost scale invariant with spectral index, ng and amplitude As- At this stage we will ignore the 
possibility of a tensor component to the fluctuations and also the possibility of early reionization, and for preciseness 
we shall assume that the temperature of the CMB is Tcmb = 2.726if , the fraction of primordial ''He is Kjjc = 0.24 
and the number relativistic degrees of freedom is — 3.04. For rest of the paper we denote 5a{trcc) — 5a, unless 
otherwise specified. 

II. VARYING a AND THE CMB 

To make a quantitative analysis of the effects of changing a on the CMB anisotropics, we to modify the linear 
Einstein-Boltzmann solver CMB FAST Here we follow the treatments outlined in Hannestad ||l^ and Kaplinhat 



et al \ VI\ and confirm their results. Changing a modifies the strength of the electromagnetic interaction and therefore 
the only effect on the creation of CMB anisotropics is via the modifications to the differential optical depth of photons 
due to Thomson scattering, 

f = XcTlcCaT , (1) 

where ionization fraction, Xe, is the fraction of the number density of free electrons to their overall number density no, 
and (Tt is the Thomson scattering cross-section. The ionization fraction, Xc, is dependent on the temperature of the 
electrons, T^, and therefore on the expansion rate of the universe a{t). Modifying a has a direct effect on the optical 
depth via the Thomson scattering cross section, ctt = Sna'^h'^ /{Sm^c^). It also indirectly effects f by modifying the 
temperature dependence of Xg. These two effects change T*, the temperature at which last scattering takes place, and 
Xeito), the residual ionization that remains after recombination, both of which influence the CMB anisotropics. 



^The astrophysical and cosmological limits that we shall discuss correspond to a limit on 5a over a particular redshift range. 
If considered in this way the OKLO constraint can be thought of as being at z ~ 0.1 
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FIG. 1. The effects of varying 5a on a standard ACDM type model as described in the text. In the top-left is the ionization 
fraction Xe and in the top-right is the visibility function -rexp[— r], both specified as functions of redshift, z. The bottom two 
curves are the angular power spectrum of CMB temperature anisotropics on the left and that of the polarization on the right. 
The solid line corresponds to 5a = 0, the dashed line to 5a = —0.05 and the dotted line to 5a = 0.05. Note that if (5q < then 
the peaks are shifted to smaller I and the amplitude of the second peak is suppressed relative to the first. 
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The change of the ionization fraction with a results from modifying the interaction between electrons and protons. 
The net recombination rate of protons and electrons into hydrogen is given by [ p9| 

t-}, - a,x,n,C - /3c (1 - ^^c) e-3Si/4Te^ (2) 

where B\ is the binding energy of the hydrogen ground state given by B\ = c?m^(? j^i^i 13.6eV for a = a(to)), and 
tte, Pe, C are constants quantifying recombination, ionization and the two photon decay in to the ground state. This 
is the net rate for recombination to and ionization from all states of the hydrogen atom. Recombination to the ground 
state can be neglected since such a process immediately creates a Lyman-a photon which reionizes another hydrogen 
atom, although one has to take into account Lyman-a photons which are redshifted out of the resonance line and also 
that the ground state can be reached by two photon decay [|9|. The recombination rate to all other excited levels 
is M 













(1) 



where A = = 2.105 x lO"^^ cm^ with the Bohr radius = h/arricC = 0.529 A and TOo = 0.511 MeV 

the electron mass. The Gaunt factor g w 0.943 is due to quantum corrections of the radiative process and is only 
weakly dependent on a. As in refs. @0 wc have ignored the effects of in change in a on this correction. The function 
<j>'{tc) = 0.5(1.735 — Into + io/6) — exp{l/tc)Ei{l/tc)/tc comes from summing up the interaction cross sections from 
all excited levels |3^. Ei is the exponential integral function and to = T^/Bi. The ionization rate is related to the 
recombination rate by detailed balance 

/3e = «c(^)'^%-^^/^% (4) 

with B2 = Bi/A the energy of the lowest lying excited, n — 2, state. The correction due to the redshift of Lyman-a 
photons and the two photon decay is given by C = (1 + KVnis) / {I + KVriis + KPcUis), with K = A^a/(87ra), 
Aa = IQnfi / {2>mcQ^ c) the wavelength of the Lyman-a photons, 2? — 8.23s^^ the net rate of the two photon decay 
with T) oc p6| , pTt and nis = (1 — Xc)'nc the number density of atoms in the Is state. 

We have incorporated these dependencies on a into CMBFAST and the results are illustrated for a simple ACDM 
model with — 0.3, h — 0.65, Vl\^h? = 0.019 and ng = 1 in Fig. |l|. On examination of the curve for Xc{z), we see 
that if the a is increasing with time {5a < 0) then the epoch of recombination is delayed, whereas if it is decreasing 
(da > 0) then recombination happens much earlier. The visibility function fexp[— r] quantifies the probability that 
a given photon observed today was last scattered at the specified redshift. Hence, one could loosely define the epoch 
of last scattering to be the maximum of the visibility function. By this definition, the ±5% shifts in a illustrated 
in Fig. |l| correspond to shifts in the epoch of last scattering by about ^100 in z, from the value of Zjoc ~ 1100 for 
Sa = 0. We also studied the effects of changing a in the process of Helium recombination and found that they are 
negligible as long as a varies in the range we discuss in this paper. 

The shift in the epoch of recombination has a number of implications for the spectrum of CMB anisotropics . 
First, the angular positions of the primary and subsequent peaks in the spectrum are determined by the physical scale 
of the sound horizon for photons at the time of last scattering. In particular, the position of the first peak is given by 
^pcak ~ T^fjo / {csV*) J where rjQ is the conformal time of the present day, 77, is that of last scattering and Cg is the sound 
speed of the photon-baryon fluid around 77*. In a model where 5a(trcc) < 0, 77* is increased while Cg is reduced by a 
smaller amount due to the larger fraction of baryons at last scattering, and rjo does not change. Hence, the first peak 
in the CMB anisotropics is moved to larger scales, or smaller €poak- Similarly, if Sa > 0, 77* and Cg are affected in the 
opposite way and the peak moves to smaller scales (larger fpeak-) A reduction of Zpoak from 240 to 200. as suggested 
by the B98 data, could be achieved by a 16% increase in Cg??*. Such effects would be degenerate in parameter space 
with modifications to Qtot, which we have ignored. 

Other effects of this shift are changes in the modulation of the peak heights by baryon drag |3^, to the photon 
diffusion damping length | |3^ , and to the time between matter domination and last scattering, which lead to subtle 
degeneracies between 6a and fib/i^ or fim/i^. The modulation of the peaks heights is determined by the relative 
density of baryons to photons at 77,, i?, = 3/3b/ i'^P-y) cx; a* oc T~^, where T*, the temperature at which recombination 
takes place, is roughly proportional to the binding energy of the electrons oc i3i oc (1 -I- Sa)'^. Hence, one might 
think that the effects of increasing the baryon density, flhh^, can be accomplished by decreasing a^ by the same 
amount. However, reducing a and delaying recombination also results in an increased diffusion photon length, which 
also could be caused by a increase in fli,h^. The degeneracy between 6a and flhh^ is, therefore, a complicated one 
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FIG. 2. The marginalized likelihoods in the Q^^h^ — h plane for the B98 (left) and B98+M99 (right) data. Shown are the 1-a 
(68% confidence) and 2-cr (95%) contours, plus a box showing the region h = 0.65 ± 0.1 and flbh'^ = 0.019 ± 0.0012 and the 
corresponding 2-cr contour. Note that the regions corresponding to the direct and CMB measurements do not overlap when 
one includes only B98 and there is only very little overlap when M99 is included as well. 

and is likely depend on the scales probed experimentally. Finally, delaying the time of recombination will alter the 
ratio of matter to radiation when the photons are last scattered, and so will have effects similar to changing fim/i^. 

III. INTEGRATED PROBABILITY DISTRIBUTIONS 

We are now in a position to calculate likelihood functions given the CMB data for Sa ^ 0. In models where a is 
constant, the B98 data, taken on their own, prefer spatially closed models with Qtot ~ 1-3 1^. In addition, many 
other parameters take similarly questionable values (e.g., to w 7.6Gyr) when the CMB data is not supplemented with 
other priors. For example, measurements of the cosmological distance ladder indicate that the Hubble constant is 
roughly between h = 0.45 — 0.85 (which we will take to be the 95% confidence range,) while the observed light element 
abundances and BBN indicate the baryon density is n^h'^ = 0.019 ±0.0024 (95% conf. level) ||l^,|lll- If 

we restrict to 

models that are spatially flat, the CMB data prefer values of these parameters in excess of that found by the direct 
measurements. 

To illustrate this point we have computed flat band power estimates for the CMB anisotropics in the range probed by 
the B98 and M99 experiments for flat ACDM models for a grid of cosmological parameters {h — 0.45 — 1.05, Ah = 0.05; 
rJb/i^ = 0.007 - 0.040, AQ.y.h'^ = 0.003; = 0.2 - 0.8, An^, = 0.1; ng = 0.8 - 1.1, Ang = 0.05.), computed the 
likelihood of the models given the data and then marginalized over the parameters ng and ri,„, assuming that (1) 
As be that measured by COBE with Gaussian errors of approximately 15% [0, (2) the B98 calibration errors are 
Gaussian, with an amplitude of 20%, and (3) the M99 calibration errors have an equivalent amplitude of 8% with 
Gaussian correlations with respect to the other measurements when included. The relative likelihood contours are 
presented in Fig. ^ for B98 alone and for it combined with M99, included also is a box giving an idea of where the direct 
limits lie. It is clear that the CMB measurements disagree with the direct measures of these cosmological parameters 
at the 2(7 level if one only takes into account B98, and this conclusion is only slightly weakened by the inclusion of 
M99. Thus, the CMB measurements appear to be in conflict with the baryon density inferred from nucleosynthesis 
and either the measurements of h or the theoretical prejudice of fitot = 1- (This is supported by the analyses of 
refs. |,§.) 

These conflicts can be resolved if one considers changing the value of a at last scattering. If we repeat the above 
analysis, there is a marked improvement in the consistency of the CMB measurements with the direct measurements 
of h and n\^h^ when Sa < 0. We illustrate this point in Fig. ^ where a is reduced by 6.5% from its present day 
value, that is 6a — —0.065. This shift brings the CMB into better agreement with the direct measurements. It should 
be noted, however, that the direct measurements themselves may also be modified by a change in the fine structure 
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TABLE L The pealis of the probability distribution functions for Sa for the various priors for the B98 data alone and for 
the combined data sets. Also tabulated are the 1- and 2-a intervals and the integrated probabihty that Sa < 0. Note that all 
the probability distributions favor Sa < 0. 



constant, which we have not attempted to model here. This issue we will discuss further in section IV 



To further quantify this, one can derive likelihoods for the value of Sa by marginalizing over all the other cosmological 
parameters, including the Hubble constant and baryon density. We have done this for a wide range of Sa using a 
slightly wider spacing than above — Ah = 0.1 and Afib^i^ ~ 0.006. In addition to the CMB data, we consider a 
number of possible prior assumptions for the parameters, particularly focusing on those involving h and S^b/i^. Results 
based on the CMB data alone, without any prior, are labelled PO. We have included two simple priors incorporating 
fairly weak constraints on the age of the universe (PI: to > 11.5 Gyr) or on the Hubble constant (P2: h = 0.65 ±0.1.) 
These tend to have similar effects, as both effectively cut off the large h region of parameter space. We have also 
considered adding weak and strong priors on the baryon density to the previously assumed Hubble constant prior, 
(P3: P2 + n^h'^ = 0.019 ± 0.006) and (P4: P2 + ilbh'^ = 0.019 ± 0.0012.) Finally, we combine the age, Hubble 
constant, and strong baryon density priors with a constraint based on the cluster baryon fraction (P5: PI -I- P4 -|- 
/b = 0.067±0.008.), where /b is the fraction of baryons to total mass deduced from X-ray observations of rich clusters. 
In each case the quoted error bars above are taken to be the 68% (1-cr) confidence level. 

The result of these marginalized distributions for Sa for PO, PI and P4 are shown in Fig. ^, and the basic properties 
of these distributions are displayed in Table I for all the priors. As can be seen, in the absence of any priors the 
data prefer a value for a at last scattering a few percent lower than its present value, but the constraint is fairly 
weak. Considering only the B99 data and adding fairly weak constraints on the age or h gives a significantly stronger 
signal, suggesting a detection of variation in a at the l-cr level. Finally, if one includes the stronger constraints that 
the baryon density is low, then the evidence for variation in a is significant at the 2-a level. Including the M99 data 
weakens these detections somewhat, but not dramatically. 

One can understand the effects of the priors in light of our earlier discussion of how the CMB anisotropics depend 
on Sa. By allowing h and Obft.^ to vary freely one can fit both the primary peak position and secondary peak heights 
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6a da 

FIG. 4. The likelihood marginalized over the cosmological parameters h, Q^, i^hh^-, ng and as a function of 5a with 
no prior (PO - solid line), the age prior (PI - dotted line) and the combined /i/strong BBN prior (P4 - dashed line). All the 
probability distributions are biased toward 5a < 0, with those for P4 implying a substantial probability that 5a < 0. 



either with (5a = or (5a ^ 0, and as we have ahea(iy seen in Fig. ^, fixing 5a — requires large values of h and 
fl\^h'^. The constraint on h comes primarily from the position of the primary peak. When one allows the time of 
recombination to change <iue to a variation in da, the high value of h can be offset by an increase in Cgfy*, making 
models with a lower values of h and i5a < equally likely. If one further assumes a prior such as PI or P2 which 
penalizes high h, the models with Sa < become favored over those with Sa — 0. A similar line of argument can 
be applied to flhh^- When 6a — a, high value of ilhh^{^ 0.031) is favored, but when one allows 6a to vary, the 
approximate degeneracy between 6a and flhh^ makes models with a lower value of flbh^ and 6a < equally likely. 
Clearly, the inclusion of a constraint which penalizes high Q,hh^ will favor models with 6a < 0. The remarkable aspect 
of this is that the inclusion of a single parameter, 6a, has the effect of improving the fit to the observations through 
two physically very distinct mechanisms. 

The conclusions are substantiated by examining the best fit models listed in Tables II and III, and plotted compared 
to the B98 data in Fig.|^ and to the combined dataset in Fig. |^. Except for PO in the case of B98 alone, and PO and 
PI in the case of B98 and M99, the reduced of the best fits are decreased by allowing a to vary, but the fits are 
only significantly improved when one assumes a strong prior for the baryon density (P3-P5). When considering just 
the B98 data, the appears to be somewhat smaller than the number of degrees of freedom, so that the reduced 
are significantly less than one, which in turn suggests that the error bars of the B98 data may be overestimated. The 
reduced no longer appear low when the M99 data are included in the analysis. 



IV. DISCUSSION 



In the previous section we presented evidence for a variation in a using recent CMB observations. If confirmed by 
subsequent observations this would be a truly remarkable result. In this section we will discuss the various aspects of 
our analysis, focusing on the potential uncertainties. 

First, we should make some comments on the details of our statistical procedure and the models which we have 
probed. We make the approximation that the CMB data points are statistically independent and Gaussianly dis- 
tributed, with window functions given by flat band powers in £. This is also an assumption in the analysis of rcf. [Q, 
but not in those of the BOOMERanG and MAXIMA collaborations |^,|| where the exact experimental window func- 
tions were used. Since the analysis presented here for 6a ^ agrees qualitatively with those other analyses, we believe 
these approximations should be sufficient for our purposes. The B98 data has more points (12 versus 10) and smaller 
errors than the M99 data, and generally provides stronger constraints. While both data sets suggest the secondary 
peak is suppressed relative to the first Doppler peak, the M99 data shows no evidence for a left-ward shift of the peak. 
Thus, including it tends to weaken the evidence for a time varying a. 
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1 1 

FIG. 5. The angular power spectra of temperature anisotropies and polarization for the best fit models to the B98 data from 
Table II. Included are PO with 5a ^ (Long dashed-dot line) and PI (Solid), P2 (dotted), P3/P4 (short dashed) P5 (short-dash 
dotted) all with 5a ^ 0, along with the B98 data. The best fit normalization changes for different models, therefore each curve 
has been appropriately normalized so that they can all be plotted with the data on the same graph. 
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TABLE II. The best fit models for the B98 data only with the various priors, with and without allowing variations in a. 
The number of degrees of freedom (number of data points minus number of theoretical parameters) for the fits are roughly 
between 7 and 10 when 5a — 0, depending on the number of constraints. (This is reduced by 1 for the 5a models.) Models 
with stronger priors included in the likelihood effectively have more data and thus more degrees of freedom. Note that the fits 
are substantially improved by allowing 5a ^ when there is a prior assumption on Q^yh^ and that the reduced a-r^ generally 
less than one. i?B98 is the ratio of the COBE and B98 normalizations for the Ce's. 



8 




1 



FIG. 6. The angular power spectra of temperature anisotropics for the best fit models to the B98 and M99 data from Table 
III. Included are PO with Sa — (dotted line) and P4 with Sa = (dashed line) with Sa ^ (solid line). With the strong 
BBN constraint on the baryon density, the 6a case provides a much better fit than the Sa — case. As for Fig. ^ we have 
normalized so that all the models can be presented on the same plot. In doing this we have fixed the calibration of the B98 
data relative to the M99 calibration to be 15% higher than its nominal value. 
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h 


fib ft'' 




ns 


Rb98 


-Rm99 


x' 


PO 


0.0 


0.75 


0.025 


0.3 


0.925 


0.92 


1.06 


16.43 




-0.025 


0.65 


0.025 


0.4 


0.925 


0.92 


1.06 


15.82 


PI 


0.0 


0.75 


0.025 


0.3 


0.925 


0.92 


1.06 


16.43 




-0.025 


0.65 


0.025 


0.4 


0.925 


0.92 


1.06 


15.82 


P2 


0.0 


0.65 


0.025 


0.5 


0.925 


0.98 


1.09 


17.13 




-0.025 


0.65 


0.025 


0.4 


0.925 


0.92 


1.06 


15.82 


P3 


0.0 


0.65 


0.025 


0.5 


0.925 


0.98 


1.09 


18.13 




-0.025 


0.65 


0.025 


0.4 


0.925 


0.92 


1.06 


16.82 


P4 


0.0 


0.75 


0.019 


0.2 


0.850 


1.00 


1.15 


21.89 




-0.065 


0.65 


0.019 


0.2 


0.900 


0.82 


0.97 


18.08 


P5 


0.0 


0.65 


0.019 


0.3 


0.850 


1.00 


1.15 


24.43 




-0.055 


0.55 


0.019 


0.4 


0.900 


0.90 


1.03 


18.96 



TABLE III. The best fit models for the B98 and M89 data with the various priors, with Sa — and Sa 0. The number of 
degrees of freedom for the fits are in the range 16-19 for Sa — 0. (Again, it is reduced by 1 when Sa ^ 0). Rmga is the ratio of 
the COBE and M99 normalizations for the Cts. 
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We include in our calculations the uncertainties in the absolute calibrations of the data sets, which is necessary in 
order for them to be consistent with each other. The B98 data was normalized by the CMB dipole, which is subject 
to large systematic errors, and their quoted calibration error is 20% in the power. The M99 experiment was also able 
calibrate off of Jupiter, and has only an 8% error. Best fit models allowing both of these to vary seem to prefer an 
increase of roughly 15% in the relative B98/M99 power calibration (see, for example. Table III). 

We should also note that the range of models we use in our analysis does not include many cosmological scenarios 
that are often considered. We have excluded the possibility of tensor fluctuations with a spectral index nx and 
amplitude At, a hot dark matter component ili,, and also that of early reionization often quantified by tr, the optical 
depth to reionization. These were included in ref. Q] and were found to have little or no bearing on the preferred 
values of fli^h^ and h since these parameters are effectively orthogonal given the present data. In fact, in ref. Q it was 
suggested that there is a degeneracy between ns and tr for the angular scales probed in B98; our unusually low values 
of ns, therefore, take into account the possibility of reionization at moderate redshift with ng = 1 being compatible 
with the data. 

We have investigated this by constructing the Fisher matrix which quantifies the effects of changing parameters has 
on the measured band powers Bi, 

where pa are the parameters (fim^^, ^^b^^, 5a, ns, tr, ^5) and dj is the data covariance matrix, assumed diagonal 
except for the calibration uncertainties. We find the strongest degeneracy of 5a to be with ilm^^, but there are also 
significant overlaps with h and Vl-i^h?", all of which are consistent with the simple theoretical arguments in Section ||. 
Our inferred errors on 5a from Figure 4 are quite consistent with those expected by computing inverse Fisher matrix 
for the best fit model with PO. The matrix is also largely block diagonal as one might have expected, with ns, tr, ^5 
being largely orthogonal to the other variables, but with significant overlap amongst themselves, confirming that ns 
and tr are degenerate given the present data. 

Knowledge of the Fisher matrix allows us to investigate the impact future CMB measurements might have on 
further constraining parameters. If the error bars of the B98 and M99 experiments are reduced by a factor of two, the 
errors on a (and indeed on most parameters) are reduced by a comparable factor. Hence, improved sensitivity with 
the same angular coverage is an important goal. We have also considered the impact of a hypothetical measurement of 
the third peak, centered at ^ = 800, as well as a detection of the first polarization peak centered at £ — 350, assuming 
10% errors on a flat band power measurement in each case. Both, particularly the polarization measurement, help to 
reduce uncertainties in f2m^^ and flhh'^, but disappointingly neither do very well in reducing the uncertainties in 5a. 
It should be noted that a weakness of this band power based approach is that the conclusions will depend somewhat 
on how the data are binned, particularly if the models vary greatly across the bins. 

Having argued that our statistical procedure and the models which we have probed provide a robust detection of 
a variation in a given prior assumptions from direct measurements of h and fl^h^ in a flat universe, we now turn 
to the more difficult issue of the effect of such a variation on these direct measurements. This pertains primarily to 
those associated with BBN since the measurements of h, to and ft, are made at very low redshifts and therefore will 
be relatively insensitive to these changes. 

Since electromagnetic effects are ubiquitous in BBN, it is clear that there must exist a constraint on 5a{tnuc) from 
the consistency of BBN with light element abundances. There are two approaches to this problem documented in the 
literature. The first is to use only the observations of ''He, which arc thought to be the most reliable. One 

can make a very simple estimate of the primordial '*He abundance, in terms of m„/mp, the neutron to proton mass 
ratio. However, expressing this in terms of a cannot be done in a model independent way because it involves a subtle 
interplay between electromagnetic, weak and strong interaction effects which have not been understood completely 
within QCD. Therefore, tight limits on 5a{tnuc) computed in this way should be treated with some caution. 

A more reliable alternative |^ is to make a detailed analysis of how changes in a can effect all the light abundances 
and derive a constraint from demanding consistency with their observed values. Although this involves a number of 
complicated nuclear reaction rates, it turns out that a model independent constraint might be possible at the level 
of around a few percent. In fact, it was suggested in ref. that |5a(inuc)| ~ 0.02 could help BBN fit the observed 
light element abundances better. However, even this value should be treated with some caution since there are further 
uncertainties which might could modify it by as much as a factor of two [^5| . 

Furthermore, the implied value of fl^^h^ is likely to change as a function of 5a(tnuc)- For example, one might expect 
that m 

^bh^ [fo(t„uc)] = [1 + Afo(W)] ^bh^ [0] (6) 
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for small values of Sa{tnuc) where the coefBcient is expected to be of order A ~ 0(1). The amplitude and sign of the 
proportionality constant A will clearly have some influence on our conclusions. In particular, if decreasing a increases 
the inferred baryon density {A < 0), then it may be possible to fit the data with a smaller change in a. However, if 
the opposite is true {A > 0), then it may prove a better fit if the fine structure constant is larger at last scattering, 
contrary to our present results. 

To make contact with the earlier discussion, one needs to relate (5a(tnuc) and (5a (tree) which requires a model for 
how the variation in a is realized. If a is increasing with time as we have suggested here, then it might be sensible to 
assume that it has done so monotonicallyQ and, therefore, Sa(tnuc) < Sa{t,:ec)- Given the uncertainties, the constraint 
from BBN on Sa{tnuc) would appear to be consistent with this relation given the fairly small values of (5a(ircc) required 
for a good fit to the CMB data. However, this certainly motivates a critical appraisal of the exact constraint on Sa 
from BBN. Including the effect of changing a on the BBN measurements would require knowledge of the parameter A 
and we have not attempted to incorporate this into our analysis. This issue should be revisited in future work when 
the a dependence of the BBN constraints are better understood. 

Finally, we should mention that realistic models in which a varies may contain one or more light scalar fields which 
mediate the precise variation. Clearly, if such a field exists it should be included in the calculation of the CMB 
anisotropics in the Boltzmann hierarchy of CMBFAST, either explicitly or as a extra relativistic degree of freedom. 
This would allow a subsequent analysis to include effects of the time variation of a, rather than just a change between 
the time of recombination and the present day. Such a field could have significant energy density during the epochs 
important for structure formation (after the time of radiation-matter equality) and it may even be possible for such 
a field to act as a quintessence field , removing the need for A. 

V. CONCLUSIONS 

The most recent CMB data provide strong evidence yet that the universe is, at least approximately, spatially fiat. 
The B98 data, however, is not entirely consistent with spatial flatness and direct measurements of other cosniological 
parameters. The situation is only slightly improved when the M99 data is included. However, if the fine structure 
constant was a few percent smaller when the photons were last scattered, then a model can be found which is consistent 
with all observations. It is clear that a change in the a is not the only possible explanation for such observations, 
and the evidence we present here could equally well be thought of as favoring other delayed recombination models, 
for example, that presented in ref. p^ . 

The evidence for a time variation in the fine structure constant is significant when a tight prior is assumed for 
the baryon density. However, the baryon density inferred from measurements of primordial abundances depends on 
nuclear physics processes at times long before last scattering at the epoch of BBN. We have argued that the values of 
6a{t,:cc) ^ —0.05 that we have deduced are consistent with BBN given the uncertainties assuming that the variation 
in a is monotonic, but that inclusion of the effect of varying a on the inferred value of fli^h^ for a given Deuterium 
abundance has been ignored, mainly due to lack of quantitative information. 

Stronger conclusions must, of course, wait for better data, such as might come from the satellite experiments MAP 
and PLANCK. In particular, these will be able to confirm whether the inconsistencies of flat models with direct 
measurements and the CMB data (such as a slight shift of features to larger scales) are real. In addition, these 
experiments should be able to break the degeneracy between a changing a and flhh^, so that a change in a can be 
tested independently of what occurred at nucleosynthesis. It is clear from Fig. |^ that the best flt models can differ 
greatly for £ > 600 in temperature anisotropics and the polarization; any experiment which probes the CMB in these 
areas will be useful in breaking these degeneracies, although the results of our Fisher matrix analysis suggest they 
will require considerable sensitivity. 

Having presented our case for a few percent variation in a at around z « 1000, it is interesting to compare to the other 
claimed detection of a change in a around z « 1, and speculate as to an explanation. Naively, 6a{z « 1000) ~ —0.01 
and Sa{z « 1) ^ — lO^'^ might suggest that Sa oc z, but this would be incompatible with BBN if extrapolated back 
to z « 10^ — 10^". In order to have a chance of being consistent with BBN, such a scenario would require that the 
variation in a terminated at some point shortly before recombination, that is, Sa{t^uc) ~ <5a(troc)- At first sight 
coupling a non-minimally to gravity via the Ricci scalar or the trace of the energy-momentum tensor as suggested in 



^In fact models have been suggested in which a oscillates [g^, although this would appear at this stage to be somewhat ad 
hoc. 
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ref. ||l8| , both of which are zero in the radiation-era and non-zero in the matter-era, might seem an attractive solution 
since the variation in a would begin at the onset of matter domination. Clearly, such ideas could have a profound 
impact on understanding of Grand Unification and this particular interpretation of the most recent observations 
presented here opens up a wide range of interesting possibilities for future research in this area. 
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